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ABSTRACT 

The microscopic approach quantum dissipation process presented by Yu and Sun [Phys. 
Rev., A49(1994)592, A51(1995)1845] is developed to analyze the wave function structure 
of dynamic evolution of a typical dissipative system, a single mode boson soaked in a 
bath of many bosons. In this paper, the wave function of total system is explicitly 
obtained as a product of two components of the system and the bath in the coherent 
state representation. It not only describes the influence of the bath on the variable of the 
system through the Brownian motion, but also manifests the back- action of the system 
on the bath and the effects of the mutual interaction among the bosons of the bath. Due 
to the back-action, the total wave function can only be partially factorizable even for the 
Brownian motion can be ignored in certain senses, such as the cases with weak coupling 
and large detuning. 

PACS numbers: 03.65, 32.80, 42.50 
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1. Introduction: The typical microscopic treatment for the model of quantum dissi- 
pation is to consider a harmònic oscillator interacting with a many-oscillator bath (e.g., 
in [1-6]). This simple and rather conventional model has been used to completely clarify 
the relation between two different approaches for quantum dissipation process frequently 
appearing in the literature, i.e., the system plus bath model and the time-dependent ef- 
fective Hamiltonian by Kanai and Calderora [7,8], since Yu and one (Sun) of the authors 
wrote down the total wave function explicitly in a form of direct product of the bath com- 
ponent and the system component [1,2]. In the discussion, because the mixed variables 
were chosen to describe the system and the bath, the wave function only manifested the 
influence of the bath on the system through the Brownian broadening of the width of 
the wave function for the system, but the back-action of the system on the bath was not 
discussed. If there indeed exists the back-action of the system on the bath, it is reasonable 
to expect that, for the individual partides constituting the bath, the mutual couplings 
among them can indirectly appear in a second order through coupling the system as an 
intermediate process. 

Another question is the relation between quantum and classical systems. In many real 
situations, the classical or macroscopic states can be represented by coherent states in the 
quantum òptics and the macroscopic quantum mechanics. Therefore it is significant to 
study how the system with an initial coherent state evolves if it really has a macroscopic 
or classical meanings and to test if it can move in the classical orbits. 

In this paper, the back-action and mutual couplings with a simple model is studied 
and manifested in the coherent representation. In the presence of both the back-action 
and the indirect mutual coupling, we also consider the meaning of the wave function of the 
dissipative system governed by the effective Hamiltonian, which is also to be determined 
in this paper. A very interesting result is that a factorizable evolution represented by a 
product of two coherent states is obtained with one component representing the system, 
which evolves according to a classical orbit and possess a classical meaning. 

2. Structure of Wave Function: This paper will consider the problem with the simple 
model consisting of a single mode boson and a bath of many bosons and its Hamiltonian 
is written as 

H = hiüb + b + ^2 hujja^aj + h + h.c , (1) 
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where £j = |^|e lcrj 's are the complex coupling constants and b + , b, and ctj are the 
bosonic creation and annihilation operators for the system and the bath respectively. This 
model can be regarded as a rotating wave approximation of the original oscillator model 
with the linear coupling J2j ÍjQ x j ~ Hj[Íj^ +a j + íj^ +a \ + h.c] of the system coordinate q 
to the bath variables Xj. 

To obtain the explicit expression for the wave function of the total system formed by 
the system plus the bath, we invoke the well-known solutions [9] 

b(t)=u(t)b(0) + Y,vAt)aj(0), (2) 
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aj (t) = e'^aM + Uj (t)b(0) + £ v js (t)a s (0) (3) 



of the corresponding Heisenberg equation, where the coeficients are 



■u(t) = exp 



~t-i(u + Au)t 



(4) 



/ . >1- exp \i (u)j — üj — Acj) t — jt/2] . . 

i>,(t) = exp (-iw,-í) HL v A J / . - 7 J , 5 



-i^í ex P [* (^j - ^ - * - 7*/ 2 ] - 1 



uAt) = -Ce-^ y 3 ~> ^ -, (6) 

v (t) [ l-exp^(^-^-Ac)t- 7 í/2] 1 

3sK ' üo + Aüü-üú s -í^/2 \ üü + Au-oüj-i^/2 J i; 

where 

{ exp[t(üjj— ^ s )t] . / 

7 ^ . (8) 
í j = s 

Here, Lamb shift Auj can be absorbed into oo to get the renormalized frequency Q = uj+Auj 
and the damping constant 7 is determined by the coupling ^and the spectrum density 
piojj) of the bath. If we chose an appropriate p(u)j) similar to the Ohmic one by Caldeira 
and Leggett [3], the above solution is exact [1,2], otherwise, it can be obtained by the 
Wigner-Weisskopf approach [9]. 

Now, we present a new method in the coherent state representation to calculate the 
evolution of the wave function in the Schròdinger picture from the explicit expressions of 
the canonical operators in the Heisenberg picture. It is different from that in our previous 
works [1,2], but quite direct and effective. If the initial state of the total system is a direct 
product |^(0)) = |0) <E> üj \4>j) an d U(t) is the evolution operator of the total system, 
the wave function |^(í)) = U(t) |^(0)) at time t can be defined by its coherent state 
representation 

*(A,{A,},í) = (A,{A,} | *(*)> = (*(0)|^(í) t |A,{A i }>*. (9) 
Here, we used the overcomplete basis 

|A, {Aj}) = |A) ® n \ X i) = N ^ M) ex P ^(°) + E V}(0) j |0> (10) 

constructed by the coherent states |A) and \Xj) for the annihilation operators 6(0) and aj(0) 
respectively. Here, the normalization constant N(X, {Xj}) = exp (— ||A| 2 — J2j ||Aj| 2 ) • 
Then, we turn to obtain an explicit expressions of U(t)\X, {Xj}) by considering the role of 
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the evolution matrix U(t) in the Heisenberg picture. In fact, since U(tyO(0)U(t) = 0(t) 
and U(t)\0) = |0) for a operator O, it is easy to obtain 

U(t)1\\,{\ j }) = N(\,{\ j })e X p U^ + E^W)) 1°) 



= N (A, {A,}) exp a(t)b\0) + £ A/3(í)a](0) |0) = |a(í)> ® ü !&(*)> (H) 



where 



/ 3 j (t) = e^ t \ j + v j (ty\+ y: vsAtyXs 

Finally, we obtain a formal factorized wave function for the total system 



(12) 
(13) 



* (A, {A,},í) = ( u(ty\ + ^ j (ty\ j ] ( e ^'% + + E ««WA. ] -(14) 



3. Partial Factorization and Effective Hamiltonian: The above wave function 
(14) is not completely factorizable because of the entanglements of the variables A and 
Xj, which are implied by the term, J2j u j(t)*^j, of the bath variables Xj modifying the 
system variable A and the term, Vj(t)*X, of the system modifying the bath one. This 
former represents the bath fluctuation due to the Brownian motion while the later the 
back-action of the system on the bath. In fact the term, Yï,j u j(t)*Xj, is caused by the 
bath fluctuation operator B(t) = Y,j v j(t) a j(fy m t ne system operator b(t), which has a 
zero thermal average, but a non-zero correlation 



<B(tyB(t') >=£4|£. 



/,(M') 



Jl 7 2 + 4 (uj - u - Au) 2 



exp 



Tiojj s 



- 1 



(15) 



The term, J2j=és v s,j{tyX s , shows the mutual interactions among the bosons of the bath 
through the system. Mathematically, if the coupling is weak with the small £j, the mu- 
tual interactions are the second order as shown by v 8 j(t)* oc £*£ s . Notice that the main 
difference between the present result and that in refs. [1,2] is the back-action and the 
mutual interactions. 

When the fluctuation can be ignored for certain cases, e.g., A is very large in the 
initial state and the coupling is weak enough, the entanglement disappears so that the 
wave function becomes a product 



N 



* (a, {a,}, t) » ( U (tyx) n 4>j (e^%- + vjityx) . 

3=1 



(16) 
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In this case, all the influences of the bath on the system are represented by the damping 
constant 7 and then the wave function is partially factorizable due to the term Vj(t)*X. It 
is not diflicult to prové that the system component (f>(u(t)*X) is governed by an effective 
Hamiltonian, which is also equivalent to the Calderora-Kanai Hamiltonian. 

To prové it, we need to return into the Heisenberg picture by dropping of the bath 
operators %(0) in b(t), namely, b(t) is replaced by b(t) = u(t)b(0). However, b(t) + and 
b(t) are not the canonical operators since [b(t),b(t) + ] = e~ 7Í . However, the Bogoliubov 
transformation gives the general canonical operators 

A(t) = ab(t) + 0b{t) + (17) 

satisfying [A(t),A(t) + ] = 1, where 

\a\ 2 - \(3\ 2 = exp (7Í) . (18) 

To given the correct Heisenberg equations for operators A(t) and A(t) + , the effective 
Hamiltonian is determined by the definition (17) as time-dependent 

H eff = ih exp (-7Í) (aa* - p(3*) A(t) + A(t) + \- {pa - de/3) A{t) + A(t) + 



-(a*(3* - (3*a*) A{t)A{t) 



(19) 



where d = à — ( / y/2+iCu)a and f3 = $ — ( / y/2—iCu)(3. Notice that the number aa*— (3(3* = 5 
should be a pure imaginary number, i.e. 5* = —5. It is not exòtic that the effective 
Hamiltonian is not unique because there is only one constraint (18). Its different forms 
correspond to different realizations of the canonical variables. 

For instance, a specific solution of Eq. (18), a = exp (7/2 + itp)t, (3 = 0, ip = ü — 
^j 2 /4 + oj 2 , gives H ef f = hexp(-ft)QA(t) + A(t) with Q = yV/4 + Cj 2 . By formally 



introducing the canonical coordinate Q = Jh/2MVt [A(t) + A(t) + ] and momentum P 



— Mílh/2 [A(t) — A(t) + ] with the varying mass M = mexp(7í), this special effective 
Hamiltonian is just of the form by Calderora and Kanai. This result can be also given 
in purely quantized version in Schròdinger picture by the direct calculation of matrix 
elements, (a| H eff \(3) = (a\ih (dU(t)/dt) U(t) + \(3). 

4. Motion of The Center of Wave Packet with Quantum Fluctuation: Let us 

consider the physical significance of the above wave function and its entanglements in 
details. In the representation of coordinate-momentum, 



g (6 + 6t) , ^-^(6-6*), 



(20) 



Xj 



- {cij + a]) , Pj = {a.j - afj , (21) 
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coherent states |A) and |Aj) are understood as the Gaussians of widths ^Jh/2uj centered 

in q = ^Jh/2u (A + A*) and x j0 = ^Jh/2uj (\j + A*) respectively. If the initial state of 
the total system is a direct product of such Gaussians, 



|*(0)> 




Qo 



n 



A, 



x 



jo > 



(22) 



the wave function at time t 



i*(*)> = iA(í)>®m^(o> 




x 



n 



defines the position evolution of center of the Gaussian wave packet 



(23) 



qc(t) = ^—(\(t) + \*(t)) 

qo exp (-i 7 í) cos (ú>) t + £ J^ {t) 

Qj(t) = exp ^— 2 7^ 2 S * n ^ a ^ * ^ ~~ ^ COS ^ ^ * 



7 . 



Sin (Wj + Gj) t + (cJj — Lü) COS (Wj + CTj) t 



(24) 



(25) 



It is known from Eq. (16) that the center of the wave packet moves along the classical 
trajectory of a damping harmònic oscillator, which is described by first term in Eq. (16) 
and perturbed by the initial displacements Xjo of the bath oscillators shown in the second 
term in Eq. (16). This fluctuation effect is just the explicit manifestation of the Brownian 
motion. For a very large initial displacement, qo, the weak coupling with small £j or a sharp 
spectral distribution of the bath with a large detuning from the renormalized frequency 
üü, one can ignore this fluctuation effect. 

Finally, we consider the back-action of the system on the bath and the mutual inter- 
action among the bosons of the bath through the motion law of the center of Gaussian 
for each boson of the bath 



*M = ^ (\j(t) + a*(í)) 



Xjocoscüjt + qo x j—Re (vj(t)) + ^ Re (v Sjj (t)) J—x s0 . 

-Ui) 



LÜ-; 



(26) 
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Notice that the back-action qo^Ju/ u>jRe(vj (i) ) is proportional to both the initial displace- 

ment qo times a Lorentz function oc [7 2 /4 + (uj — Cu) 2 ] \ Thus it can not be neglected for 
large qo or in the quasi-resonance case that the bath spectral distribution p(ujj) peaked in 

the renormalized frequency ü of the system. The last term J2 s (^j) ^ e ( v s,j (t) ) \J & s /ujX s o is 
of second order and explicitly reflects the mutual coupling among the bosons of the bath, 
which can be neglected in the first order approximation. 

5. Summary and Discussions: In summary, we first mention that the Langevin ap- 
proach is an Standard treatment for quantum dissipation process in the present model 
Hamiltonian [4], but it hardly concerns the structure of wave function that is essentially 
important in zero temperature; the Markoff approximation is also a quite effective den- 
sity operator method, but it only considers very a few intuitional picture based on the 
classical correspondence and its dissipation-fluctuation relation directly. This paper take 
both two aspects of this problem into account and thus gives a direct and clarified picture 
for quantum dissipation process with the most simple model. Our discussion not only 
concerns the necessary details in the dynamics of quantum dissipation, but also reveals 
the roles of the back-action of bath and the mutual coupling among the bosons of bath. 
Notice that the Langevin approach is based on a stochastic equation 



b(t) 



1 

-- - i (u a + Au) 



b{t) + F{t) (27) 



with the stochastic force F(t) = — i J2j Íj a j{^) e ~ luJjt - Its explicit solution is the starting 
point finding the partially factorizable wave function of the total system in this paper. In 
this sense, our study is a generalization of quantum Langevin theory. As for the Markoff 
theory, our explicit solution somewhat can be regarded as the zero-temperature result of 
the density matrix approach, but it deal with the dissipation - fluctuation relation in the 
framework of wave function. It is straightforward to extend present discussion to the finite 
temperature case. Another emphasis of this paper is the partially factorizable structure 
of the evolution function, which clarifies the new meaning of the wave function governed 
by the time-dependent effective Hamiltonian by Calderora and Kanai in the presence of 
the back-action and the indirect mutual interactions. Notice that this meaning in our 
previous works [1, 2] was only understood in the sense of the complete factorization of 
wave function. 

To conclude this paper we point out that the method used both in this paper and 
the previous works [1,2] is very limited to a linear system we considered before, such 
as the harmònic oscillator, the inverse harmònic oscillator (the harmònic oscillator with 
image frequency, u — > iu) and a linear potential for a constant force. This is because, 
only for these systems, the solutions of canonical Heisenberg operators are the linear 
combinations of the system variable and the bath variables and thus the wave functions 
are factorizable in the appropriate representations. There still exists the difliculties in 
principle to generalize the idea and method of this paper to the nonlinear cases except 
those systems that can be reasonably linearized. 
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